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and subgrid flux for filtering schemes used in the turbulence modelling method of Large 



We continue a previous work in which a priori estimates were derived on subgrid stress 
1 subgrid flux for filtering schemes used in the turbulence modelling method of Large- 
Eddy Simulation (LES). The estimates were derived there as rigorous consequences of the 



o 

, exact subgrid stress formulae from Navier-Stokes equations under the conditions assumed 



for velocity fields in the Parisi-Frisch "multifractal model." It was also shown that these 
assumptions are realistic in an extended inertial range. Therefore the estimates must be 
obeyed by any faithful subgrid model and we use them here to evaluate some popular 
models of the subgrid stress (Smagorinsky, Bardina, etc.) We also examine the effects of 
the choice of filter function on the magnitudes of subgrid stress and transfer. Under mild 
assumptions on the filter these quantities are determined by local-in-wavenumber, inertial- 
range interactions and can be modelled in a universal way. However, one common choice 
of filter — the sharp cutoff filter in Fourier space — does not satisfy the modest required 
conditions and we show that the associated flux, including "backscatter" effects, may be 
spuriously dominated by nonlocal-in-wavenumber, convective processes of a non-universal 
type. 
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1 Introduction 

In the previous part of this work JjJ] (hereafter denoted I) we derived a priori estimates on 
subgrid stress and subgrid energy flux as they appear in the filtering approach employed in 
large-eddy simulation (LES). The paper I was devoted to more purely theoretical issues. The 
relationship of the estimates to Kraichnan's "refined similarity hypothesis" for local energy 
flux [Q] was also discussed, in our earlier work H. It was shown that the estimates arise from 
the fact that LES flux with a "good" filter, concretely specified by modest conditions, gets its 
dominant contribution in the inertial range from local wavevector triads. Thus the traditional 
picture of the local energy cascade was verified and an alternative scenario of "local transfer by 
nonlocal triads" suggested by DNS studies [Q, [f| was shown to be an artefact of the measure 
of energy flux employed in those works based upon a sharp Fourier filter. The sharp cutoff 
filter is the most common example of a "bad" filter, which violates the modest conditions we 
imposed. In that case, we showed that the associated energy flux may indeed be dominated 
by distant wavevector triads, whose physical origin is transfer of energy back and forth across 
the sharp spectral boundary at inertial-range wavenumber k due to rapid convection of eddies 
of size ~ 1/k by energy-range eddies. However, it was pointed out that energy flux defined 
by a "good" filter is a scale-average of the sharp spectral flux. This leads to cancellation of 
the nonlocal contribution because of detailed energy conservation instantaneously in each flow 
realization, as already established in our earlier work for a special example ||. The LES flux 
with a "good" filter was therefore shown to give an appropriate combination of wavevector 
triads in which convective effects to transfer cancel, concretely realizing independent proposals 
of Waleffe and Zhou || that such a combination should exist. All of the previous results 
were derived based upon the Parisi-Frisch "multifractal model" of inertial-range velocity fields, 
which assumes local Holder exponents h < 1 distributed throughout the flow as the molecular 
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viscosity goes to zero |lO|, |TT| . It was shown that such a hypothesis is necessary, with further 
h < 1/3 on some space set, in order to maintain constant mean energy flux. Additional rigorous 
results of [jl2| were reviewed which establish that the local Holder continuity must hold almost 
surely in every realization under very widely accepted assumptions, such as power-law decay of 
the mean energy spectrum. 

In this paper we will discuss the implications of our estimates in I for practical LES work. 
Let us just briefly recall the motivation of this approach |l3|, 14, [T^]. If a filtering function 
Gt(r) = l~ d Gz{r/l) is convoluted with the velocity field, as 

v/(r) = J d d r' G^r-r')v(r'), (1) 

then this provides a convenient definition of the "large-scale velocity" or "resolved field" v^. The 
"small-scale velocity" or "subgrid field" can be simply defined as the complementary component: 

vJ(r) = v(r)-Vj(r). (2) 
When the filtering operation is applied to the Navier-Stokes system an equation is obtained for 

d t vg + V • (v£V> + n) = -S7pn + vq A vg, (3) 

in which is the filtered pressure field (required to maintain V • = 0) and ti is the subscale 
turbulent stress. It is given as an explicit quadratic function of the velocity 

r^ = T £ (v,v) (4) 

with 



T>(u, v) = (uv) ; - u^. (5) 

Because of this last term the Eq.(||) is not a closed equation for v^. The main effort of LES 
practitioners is to devise suitable model representations T£ = T™ odol (v^, v^) as functions of 
in order to close these equations, allowing them to be solved numerically. (It is actually 
only the trace-free or deviatoric part, TV= Te — 3 (TrT^) I which is generally modelled, while 



the diagonal part is absorbed into the pressure term.) The theoretical justification for such 
modelling lies in the presumed universality of the small-scales v' e when the grid-scale parameter 
£ is chosen as an inertial-range length. In that case, a generally applicable or universal model 
of the stresses may have some validity. The computational advantage of LES in very high 
Reynolds number flow is that the filtered Eq.(^) may be solved numerically on a computational 
grid with point-separation A ~ 1/2 rather than with a mesh which resolves scales down to the 
viscous length rj, as necessary for direct numerical simulation (DNS) of Navier-Stokes. It is 
more customary in the LES literature to denote the filtering scale as A rather than £, because 
of this fact. 

One of our main theses here is that the estimates we have established may be used to place 
some constraints on models a priori since the bounds satisfied for the exact subscale stress 
must be satisfied also for a faithful model. That is, given the true velocity field, with local 
Holder regularity in the inertial-range, the two functions T^(v, v) and T™ odel (v^, v^) should 
scale in the same way at each point r as £ — ► through that range. Although an actual 
model simulation is performed with fixed £ = A, it seems reasonable to require the previous 
correspondence for £ — > since it is based upon the ideas of universality and scale-similarity 
which underlie the whole scheme. In particular, a model which has this property will faithfully 
reflect "intermittency" effects due to the range of Holder indices h < 1 in the flow associated to 
structures either more or less singular than the K41 mean- field value h = ^. This is correct at 
least for the model stresses and true stresses obtained by filtering the true velocities, if not for 
the model stress obtained from the LES simulation itself compared with the true stress. The 
proposal is similar in nature to the so-called a priori test which was first suggested by Clark et 
al. (and sometimes called the "Clark test"). In that case, it was proposed to judge models 
by their ability to agree with the true subscale stress calculated by filtering velocity fields 
obtained as data from a (nominally) high Reynolds number DNS or well-resolved experimental 
measurements. It has also the same spirit as the "statistical a priori tests" proposed recently by 
Meneveau [17]. As for both of these proposals it is recognized that satisfaction of the a priori 



test cannot guarantee that a model will perform well for some specific task, such as accurate 
prediction of mean velocities or energies, faithful reproduction of large-scale flow structures, etc. 
On the other hand, if the a priori estimates are invalid for a specific model, then it shows that 
the model will give erroneous predictions for the true stress even if the filter is applied to true 
velocities. Hence errors will be incurred in addition to those that arise from the discrepancies 
between the true resolved velocity and that produced by the numerical LES. For that reason the 
satisfaction of the a priori estimates by the model stress is a necessary condition to accurately 
reproduce the subscale stress in a numerical LES, unless some fortuitous cancellation of errors 
can occur between the false scaling of the model and the deviations of the true and numerical 
velocities. 

A second main thesis of this work involves the effects of the filter choice in LES work. It 
was already observed in I that very distinctively different results are obtained for magnitudes of 
energy transfer in the case of "good" filters which satisfy our modest criterion and "bad" filters 
which violate it. This point has practical importance since one commonly used filter, the sharp 
Fourier cutoff filter, is "bad" according to our criterion. In that case, the convective effects 
discussed in I will be present if the cutoff filter is applied to the actual velocities. In particular, 
as we discuss further below, for a "bad" filter the distant triads produce — in the usual language 
of LES — a large backscatter of energy which is convective in origin, whereas with a "good" 
filter the backscatter is due only to local-in-wavenumber straining processes. We therefore draw 
here an important distinction between convective backscatter and strain backscatter, which does 
not seem to have been emphasized before. Since the convective backscatter has its origin in 
interactions with the large-scale, inhomogeneous eddies, it cannot be universally modelled. Any 
attempt to do so will require an elaborate alteration of the model tailored to the large-scale 
motions if it is even possible at all. This type of problem, which is predicted by our rigorous 
a priori analysis of the equations of motion, will be shown to occur in the previous attempts 
to model "convective sweeping," such as that of Chasnov Jig ]. Hence we shall argue that 
modelling of the convection backscatter is ill-advised and, practically speaking, impossible for 
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inhomogeneous flows. Even for homogeneous turbulence it leads to reduced computational 
economy of the simulation without any corresponding increase in its faithfulness or accuracy. 
For a "good" choice of filter the convective backscatter is entirely absent and hence there is no 
rationale to incorporate it in numerical LES. One should keep in mind that in most current 
LES computations the explicit filter function never appears anyway unless a "defiltering" of 
the computed velocity is attempted. Therefore, the primary practical implication of this work 
is on the choice of filter used for filtering studies of either DNS or numerical data to determine 
what statistics of the subgrid stress ought to be used for modelling purposes. 

The precise contents of this work are as follows: in Section 2 we shall derive a somewhat more 
general form of our estimates than in I and discuss their validity for specific filters commonly 
employed in LES. Thereafter we shall discuss implications of the estimates for LES and, in 
particular, investigate their a priori validity for some common stress models. In Section 3 we 
shall discuss the issue of "convective vs. strain backscatter" based upon our previous work in 
I and explain why the former must lead to non- universal effects. We then discuss how such 
behavior is seen also in analytic closures of EDQNM-type and in previous attempts based upon 
them to model "convective backscatter." Morals for LES will be drawn. Finally, in the Section 
4 we will make our conclusions. 



2 A Priori Estimates and Evaluation of Stress Models 

(2.1) Derivation of the Estimates 

We shall succinctly rederive and also generalize the estimates in I for Tg. To that end, let us 



review here the Germano decomposition of the subscale stress [19], which also forms the most 
rational basis to introduce the common subgrid stress models in the next subsection. The idea 
of the decomposition is straightforward. One substitutes v = + \\ into Tg to obtain: 

T<(v,v) = L/(v,v) + C/(v,v)+R«(v,v), (6) 
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where 

L e (v,v)=T t (y lf V e ), (7) 

is the so-called Leonard stress, 

Q(v,v) = T t (v t y t ) + T i (y' i ,v e ), (8) 

is the cross stress, and 

R £ (v,v) =T e (V e y e ), (9) 
is denoted as Reynolds stress. To be more precise, these are "modified" forms. 1 In the 



Germano decomposition, as its originator himself emphasized |l9j], each term is separately 
Galilei invariant, whereas this is not so for the traditional decomposition. 

We now make here a similar — and related — observation: each term separately in the Ger- 
mano decomposition Eq.@ obeys the a priori estimate derived in I for m alone. We recall that 
the estimates in I were based upon the local Holder condition postulated in the Parisi-Frisch 
"multifractal model" : 

|v(r + l)-v(r)| =o(f) , (14) 

with < h < 1, which is usually denoted v G C h (r). The reader should refer to I for the 
theoretical basis of this assumption. Under this hypothesis it was proved in I that 

T/(v,v)=oM (15) 



1 The traditional decomposition |2C|] was instead 

T e (v, v) = tt(y, v) + C e (v, v) + R e (v, v), (10) 

with 



L<(v, v) = (vtvt) t - ViVt, (11) 

and 



Q(v,v) = (v^), + (v^)„ (12) 

and 



R,(v,v) = (v'v').. (13) 
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at the space point r. Here we shall extend this same estimate to the separate terms L^, C^,R^. 
This follows as a consequence of two facts. 

First, there is a simple generalization of the identity of Constantin et al. pl| used in I. It 
reads 

T*(u, v) = (Au(r)Av(r)}, - (Au(r)) { (Av(r)) f , (16) 
where Aiv(r) = v(r) — v(r — 1) denotes the "backward difference," and 

</>< = J d d l G,(l)/(1), (17) 

is a "separation-average" of a function / of the displacement vector 1 with respect to the filter 
function G#. The proof of it is the same as in I. We may note also that 

vKr) = (Av(r)) (18) 

and 

Vv £ (r) = - J d d \ (VG,)(l)A lV (r), (19) 

which were proved in I. The first identity, Eq. (|i~6|) , has been independently discovered by Vreman 
et al. |22]| , who employed it to establish realizability inequalities for the stress tensor. Note 
that this identity also makes the Galilei invariance of each separate term obvious, since each is 
expressed entirely in terms of velocity differences. 

The second fact used to extend our estimates is the observation that, if v £ C h (r), then 
also 

Ajv^r) = O (f) , (20) 

and 

A 1 v' e (r) = 0{£ h ). (21) 
A sufficient condition on the filter for these results to hold is 

J dx [|G(x)| + |VG(x)|] • |x| 2 < oo, (22) 
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which was the condition for a "good" filter set forth in I. It is not necessary and may actually 
be weakened to some extent, but it is anyway very mild in practice. Eq.(|20|) follows directly 
from the identity 

Aiv,(r) = Jds G e (s) [A 1+S v(r) - A s v(r)] . (23) 

Note that derivation of the estimate Eq.([20|) here uses only the weaker condition / e?x |G(x)| < 
oo on the filter function. Since 

A,v$(r) = A lV (r) - Ajv^r), (24) 

the second estimatee Eq.(pT|) immediately follows as well. 

Now we can easily see that the same big-0 estimates hold for each of ~Le,Cg, and R.£ as 
for the entire tensor T^. Because each of the terms in the Germano decomposition just comes 
from replacing each v in Tg(v, v) by a v> or a v^, our earlier analysis of I gives the estimate 
separately to each term. Note that the same is not true of the terms in the traditional 
decomposition. This is directly related to their lack of Galilei invariance, since each term can 
have a proportionality to the amplitude of the homogeneous velocity component. The validity of 
the estimate implies that each of the terms scales as ~ v^v^, just as suggested by conventional 
heuristics. In particular, it follows under our modest filter requirement Eq.(|2"2j) that the stress 
components are truly small-scale quantities, locally-determined in scale and without direct 
dependence upon the large-scale modes. It is precisely this fact that makes possible a universal 
model of the stresses, generally applicable to any flow. On the other hand, the cross and 
Leonard stresses in the traditional decomposition scale instead as v^v^. The "subgrid stress" in 
that case is not truly a subgrid-scale quantity but depends also upon the large scales of motion. 

The estimates on the subgrid stress have another interesting implication which was pointed 
out to us by U. Piomelli. Since the Holder condition must be satisfied in the inertial-range of 
turbulent flow, then Eq.(|i~5|) gives some justification to the use of 2nd-order accurate numerical 
schemes in integrating LES models. In that case the errors are formally 0(£ 2 ), which is asymp- 
totically negligible compared to ti ~ £ 2h for small £, when h < 1. Reynolds had earlier made an 
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estimate that m ~ £ 2 and argued that a second-order scheme would entail errors of the same 
order as the stress itself and, therefore, be inadequate (U. Piomelli, private communication). 
However, Reynolds' estimate assumes the validity of a Taylor expansion to first-order for the 
velocity field, or h > 1. We have seen in I that this, or even a much weaker condition h > |, is 
inconsistent with constant mean flux of energy in the inertial range. Therefore, it is not a valid 
estimate in the inertial-range and second-order schemes should be numerically adequate. 

Let us consider the validity of the condition Eq.(p2[) for filters commonly used in practical 
LES. The Gaussian filter is a "good" example which easily satisfies these constraints. It is 
usually defined in the LES literature as 

6|x| 2 " 



G(x) = -) exp 
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(25) 



However, two other common choices do not satisfy our weak condition. The tophat filter 

1 if max,- \xj\ < h 
(rW = { 2 (26) 

otherwise 

is not differentiable. Nevertheless, although our sufficient condition is violated, all of our esti- 
mates in I and here are still valid for the tophat filter. The only place where the differentiability 
of the filter was used was in deriving the estimate 

Vv,(r) = O (V 1 " 1 ) . (27) 



This still holds for the tophat filter, since 



v«w = F n/_ 



i i t — r r 1 / 2 

dh k 



1-1/2 



v | r + h k e k + ~e^\ - v + h k e k - ^ 



(28) 



with that filter choice. The sharp cutoff filter also violates condition Eq.(|22|), this time because 
of a more serious problem of slow spatial decay. The filter is defined by its Fourier transform 



1 if maxj \kA < ir 



G(k) = { (29) 
otherwise. 
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In physical space this gives 

G{x)= 2 d f[ a -^1, (30) 
i=l 



sin irxi 
1=1 '' ' 

-J 



which goes as |x| for large x. Because of this exponent-d power-law decay, the moment 
condition is not satisfied. In fact, we will see by example in the next section that the previous 
bounds for T£ may actually be violated for this filter. 
(2.2) A Priori Evaluation of Stress Models 

We are now in a position to introduce and evaluate the common subgrid stress models 
according to the criterion discussed in the Introduction. 

The most common models are very easily described on the basis of the Germano decompo- 
sition. Note that the first Leonard stress term is already a function of the resolved field Vi and 
does not need to be modelled. The remaining terms, however, involve and must be modelled. 
The simplest, and oldest, model is the Smagorinsky model [23] which takes 



T e SM (V e ,V e ) = -2u t W e , (31) 
where u^(r) is a local "eddy- viscosity" determined by the formula 



v t ee {Csty pvj, (32) 

which arises from an energy-balance argument. Within the Germano decomposition, this may 
be thought to come from dropping the Leonard term and taking the previous model for the 
C + R stress terms. Exactly the opposite choice is made in the Bardina model [24], which drops 
the cross and Reynolds terms and simply keeps the explicit Leonard stress L: 

Tf B (v e ,V e ) =T(yi,v e ). (33) 

The most intuitive procedure in these terms is, of course, both to keep L and to model C + R 
a la Smagorinsky, which yields the mixed model: 

TV MX (v £ , Vt) ee T(v £ , V e ) - 2v e * e , (34) 
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These may be regarded as the "classical" subgrid stress models for LES. 

It is very easy to verify that these models all satisfy the a priori estimates given our prior 
discussion. In fact, since T^ D = identically, the estimates for the Bardina model follow from 
our discussion above for the Leonard stress. In the case of the Smagorinsky model we infer 



from Eq.(j27|) that 

a e = U^ 1 ) , (35) 



which gives also 

v t = (l h+l ) (36) 



from the prescription Eq.(32) for the local eddy viscosity. Therefore, it follows that 

T e SM = o(f h ), (37) 

as required. Finally, the mixed model is just the sum of the previous two, so it obeys the 
estimates as well. As a consequence, all of these models will, if evaluated for the true velocity 
fields in the inertial range, obey the same big-0 estimates as the true subgrid stress at the same 
space point. This gives one some hope that the models will be able to capture "intermittency" 
effects due to structures in the flow. 

Let us discuss some other models besides the above "standard" models. A rather recent pro- 
posal, made and implemented by Germano et al. | j25| |, attempts to improve upon the Smagorin- 
sky and mixed models by using information from the evolving flow itself to fix the parameter 
Cs, which is no longer a constant but ranges over both positive and negative real values. This 
is the so-called dynamical model. We shall not discuss all the versions of this model — of which 
there are several in the literature — or its precise derivation. In none of the published cases 
have we been able to prove the a priori estimates. However, this is not necessarily a nega- 
tive conclusion. Merely having unwieldy expressions for which it is inconvenient to establish 
a mathematical estimate is no reason to disqualify a model. Rather it must be shown that 
the estimates are definitely violated. For example, it may be possible to verify the estimates a 
posteriori. In the case of the dynamical model the estimates will be satisfied according to our 
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previous analysis if the model parameter C5, while not constant, stays between fixed bounds 
as i is decreased. In the work of Zang et al. on a lid-driven cavity flow pq ], the parameter Cs 
was observed in the midplane of the flow to range from -0.12 to 0.1 for a Smagorinsky-based 
model and from -0.018 to 0.026 for a dynamical mixed model. If these ranges of values were 
observed not to increase under successive refinement of the grid, then the estimates would be 
a posteriori verified. This, then, is an important question to consider numerically. 

So far we have not discussed any "bad" models. Some of the earlier models based upon the 
traditional Leonard decomposition of the subscale stress, e.g. the model used by Moin and Kim 
in an early LES study of turbulent channel flow |p7[| , have problems in this regard. We shall 
not discuss this in detail, since these models are no longer in common use, but we shall just 
mention that the problems with the estimates in these cases are directly connected with the 
lack of Galilei invariance discussed by Germane However, there are more subtle things that 
can go wrong. All of the estimates verified above depended upon using a "good" filter that 
satisfies a modest condition like Eq. (^2|) . Nevertheless, there are "bad" filters for which the a 
priori estimates are not satisfied, including one popular choice, the sharp spectral cutoff. We 
now discuss some of the problems this entails. 

3 The Fourier Cutoff Filter and Non-Universal Effects 

(3.1) Violation of the Stress Estimates 

We have already observed in I that the energy flux defined with the sharp spectral filter 
does not satisfy our a priori estimates because it is dominated instantaneously by nonlocal, 
convective interactions. It is not enough to observe that the condition Eq.(^) is violated, 
because it is only a sufficient condition and not a necessary one. Instead, the estimates must 
be directly verified to fail. For that purpose we used an example from || , of the form 

v(r) = 2f7(ei - e 2 ) sin [(ei + e 2 ) • r] 
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2A c r i r i ~i 

+ E {cos [2 w e! • rj + cos [(2 N+1 - l)ei • r - e 2 • rj } , (38) 

7V>1 

with U,A real constants and < h < 1. This velocity field is given by a so-called Fourier- 
Weierstrass series and it is known that v € C h (r) for every space point r in the periodic box 
where it is defined. We will just cite the result obtained in Q, without repeating the calculation, 
that 

II A = (2 2 - h UA 2 ) x 2 A ( 1 - 2fe ) (39) 

where IIa is the (global) subgrid energy flux across wavenumber 2 A defined with the sharp cutoff 
filter. Observe the very important feature that the flux is here proportional to the amplitude 
U of the lowest excited wavenumber mode. Note that this is not due to any failure of Galilei 
invariance. The spectral flux H(k) for k ^ is perfectly invariant to changes in the zero- 
wavenumber amplitude. However, its leading term may still have a direct dependence upon the 
amplitude of very small wavenumber modes for arbitrarily large k. The physics of this, as we 
discussed previously, is the back and forth transfer of energy across the spectral boundary k by 
small steps, associated to the rapid convection of small structures of scale ~ 1/k by the largest 
eddies. In the language of LES modelling there will be a large backscatter due to this effect, 
which we may call convective backscatter. 

It might appear that the example we have considered is rather contrived. However, a little 
thought shows that, so long as one keeps a sharp spectral cutoff, the same types of effects will 
occur generically in Ug . In fact, if a given velocity field had instantaneous transfer determined 
by local triads, its superposition with the above velocity field will have transfer dominated by 
the convective, distant triad interactions, which are asymptotically much larger. The physical 
interpretation also argues that such effects will be quite typical. However, we have seen that this 
type of process — while a real physical effect — is essentially spurious in terms of energy transfer, 
since it averages away over time and cancels instantaneously in sensitive measures of the flux. 
In particular, we showed in I that this contribution from distant triads cancels in LES flux with 
a "good" filter, due to detailed energy conservation. The reader is encouraged to review Section 
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3.1 of I where this cancellation of the convective contribution was demonstrated. It must be 
distinguished from the backscatter which arises from the local, straining processes whose contri- 
bution to energy flux in 3D is forward on average — because stretching of isotropic distributions 
of small-scale vortex elements by large-scale strain increases their mean energy — but which also 
has both negative and positive fluctuating values. Since these local, stretching interactions are 
the crucial ones for energy transfer, it is this second backscatter, or strain backscatter which is 
more meaningful to the dynamics. Nevertheless, the "convective backscatter" has, as we have 
seen, a much bigger absolute magnitude than the "strain backscatter" since it is nonlocally 
driven by the largest, most energetic eddies. Therefore, if an inappropriate filter is used — such 
as the sharp Fourier cutoff — then the "convective backscatter" will completely dominate and 
mask the more dynamically intrinsic "strain backscatter." 

This seems to be exactly the effect observed in numerical simulations of Piomelli et al. 
It was found in that work that with the sharp cutoff filter the fraction of space points 
instantaneously experiencing backscatter was about 50%, compared with about 30% for the 
Gaussian filter in the same flow. Furthermore, the intensity of backscatter with the sharp 
Fourier filter was as much as an order of magnitude higher than for the Gaussian. Both of 
these facts are consistent with the analyses we have made and point toward large convective 



effects in the sharp cutoff case. It was independently argued by Vreman et al. [22] that the 
non-positivity of the cutoff filter might increase the fraction of space points with negative flux. 
This argument is not so clearly valid, since the flux is not positive even with a positive filter and 
can take either sign. Our explanation is rather that the large convective contributions to the 
flux of both signs, positive and negative, are present with the cutoff filter. Since the convection 
effects will cancel on average even for the cutoff filter, one should expect that the fraction 
of points with either sign of flux will be roughly 50% — 50% and that the overall magnitude 
will be much higher. Both of these effects are observed. Notice that our explanation of this 
effect has predictive power. In fact, if the large backscatter observed in p8| is "convective" 



in origin, as we have proposed, then it ought to be absent if the same study were performed 
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for Kraichnan's "modified Navier-Stokes system" (Section 6 of p9|| .) That dynamics lacks the 

nonlocal, convective interactions which we have proposed as responsible for the effect, so that 
the results obtained there for the Gaussian and the sharp Fourier filters ought to be essentially 
indistinguishable. This is a crucial test that ought to be performed. 

Even if our interpretation is accepted as correct, it might be thought by some that there is 
nothing intrinsically wrong with the sharp Fourier cutoff. Some account must be taken of the 
difference in physics for the "good" and "bad" filter cases — exemplified by Gaussian and sharp 
Fourier, respectively — but a subgrid model with either of these filters could be constructed so 
long as it were consistent with the processes that actually occur with that filter. Since the 
"convective backscatter" actually occurs with the sharp cutoff filter, a faithful subgrid model 
should be constructed to represent it. However, we emphatically do not agree with such an 
assessment. It is crucial to understand the source of the convection dependence. If one uses 
the sharp Fourier filter, then the local flux is still represented by the formula ILj = -Vv^ : ti 
just as for any filter. However, it is not too hard to show that ||Vy«||oo = O (i 11 " 1 ^ here just 
as for the "good" filter case. (For example, it follows from the arguments that establish the 
"Littlewood-Paley-type" conditions for Holder continuity: see ||.) Therefore, the violation of 
the estimate for the energy flux with the sharp spectral cutoff is due to a violation of the estimate 
for the "subscale stress. " In other words, the nonlocal, convective type interactions must appear 
in this "subscale stress" so that it is ~ £ h rather than ~ £ 2h . The expected estimate T£ ~ v ^ v £ 
is incorrect and instead T£ ~ v^v^. 

Such anomalous behavior is explicitly verified for the example above, as we now show. Our 



between the top mode of 5a and the bottom mode of Sa+i- To simplify formulas, we introduce 
some notation, as follows: 



calculation will be done with ^ 



k chosen to lie in the interval 




U = J7(ei - e 2 ), A = Ae 3 



(40) 



ip = (ei + e 2 ) • r, 9 N)1 = 2 e x • r, 6 N>2 = (2 



,N+1 



l)ei • r — e 2 • r. 



(41) 
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Note that 

A 

v £ (r) = 2Usin<p + 2A 2~ Nh [cos 6 N<1 + cos 9 N>2 ]. (42) 

N=l 

It then follows from an elementary but slightly tedious calculation that 

L<? = -(v£vi)i = 

-2(UA + AU) • 2- Ah sin 6> A +i,i 
A 

-2AA £ 2 -( N+ ^ h [2cos(# A , 2 + e N>1 ) + cos(0 Ai2 + 9 N , 2 )] 

N=l 

-4AA • 2~ 2Ah cos A +i,i, (43) 

and 



-2(UA + AU) • 2^ A+1 ^sin# Ai2 
A 

+2AA 2-( JV+A+1 ) /l [2 cos(0 A +i,i - 9 N , 2 ) + co6(0 A +i,i - ^,1)] , (44) 

N=l 

and 



00 

+4AA ]T 2-( 2N -^ h ■ cos 99 

N=A+2 

00 

+4 A A ^ 2~ 27V/l . (45) 

Af=A+l 

As one would expect, the Reynolds stress is here 0(2~ 2Ah ), the naive result. However, the 
Leonard and cross stress terms, and C^, are not, but give contributions ~ v^v£. In fact, 
the dominant terms are ones ~ U ■ A2~ Ah , proportional to the amplitude of the lowest excited 
mode. What the estimates mean is that the "subscale stress" for the cutoff filter is not a 
true function of the small-scales alone — as one would naively believe — but there is a residual, 
implicit dependence upon the large-scales. In principle there is no difficulty with the Leonard 
term, because it does not need to be modelled and can be calculated explicitly. However, the 
cross term is a real problem. 
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This fact has severe negative implications for LES modelling with such a "bad" filter as the 
sharp Fourier one. The entire justification for the LES scheme is the universality of small-scales 
in the inertial-range, which allows a model description of the subscale stress to be attempted. 
However, if the "subscale stress" for a particular choice of filter retains a dependence upon the 
low-wavenumber, or large-eddy, degrees of freedom that are highly nonuniversal, anisotropic, 
etc. then we see absolutely no reason to believe that a simple, general model could be con- 
structed. One may ask how this should affect actual LES practice. After all, in most current 
LES computations the explicit filter function never appears anyway unless a "defiltering" of the 
computed velocity is attempted. Therefore, the primary practical implication of this work is 
on the choice of filter used for filtering studies of either DNS or numerical data to determine the 
statistics of the subgrid stress for modelling purposes. Indeed, the key implication of our analy- 
sis is that the large backscatter observed with the Fourier filter in studies such as Jjffi/ ought not 
to be incorporated into LES models, since it is just an artefact of that pathological filter choice. 
It is not even possible to model such backscatter on the basis of Kolmogorov inertial-range 
ideas, since its origin is nonlocal and is dominated by non-universal, low-wavenumber modes. 
To clarify these points, we consider now some specific models. 
(3.2) Evaluation of Models Based on the Spectral Cutoff Filter 



All of the problems discussed appear in the stochastic LES model developed by Chasnov [18], 

who used the sharp spectral filter in conjunction with a Langevin model of the EDQNM closure 

equations. From our foregoing discussion it follows that — in so far as the Langevin model is 

consistent with reality (the Navier-Stokes equations) — use of the sharp Fourier filter will lead 
2 However, this is not the case in some more sophisticated approaches, such as the "dynamical model," where 
the filtering function appears explicitly. Also, on theoretical grounds, filtering the velocity at successive time- 
steps in the course of the computation ought to be done in LES computation, with a filter length I at least as 
large as 2 x A, the grid-length of the numerical discretization. This would guarantee that the computed velocity 
field represents at every step a properly filtered object and does not develop unrealistically large energy in the 
grid-scale modes. Needless to say, this prescription is not usually adopted because available numerical power 
makes it impractical. 
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to forward (eddy viscosity) and backward (backscatter) contributions to energy transfer which 
separately contain the strong effects of convection by energy-scale eddies that only cancel in 
the net flux. In fact, Chasnov observed exactly this in his model and, indeed, claimed it as a 
virtue of his approach that he modelled "the random sweeping of small scales by large scales." 
We do not agree with his appraisal. Our work has shown that this type of effect is not present 
if a "good" filter is employed for Navier-Stokes dynamics and need not then be modelled. We 
should stress that we do not disagree with Chasnov on the issue of stochastic modelling of 
backscatter. Even with a "good" filter what we have termed "strain backscatter" will still be 
present and represent an important part of the transfer process. It is only the "connective 
backscatter" which we have proved unnecessary to model. It is clearly preferable from the 
point of computing economy to dispense with modelling such effects which are just artefacts 
of a "bad" filter choice. Most importantly, any attempt to model such effects severely reduces 
the universality of a model, since it requires explicit consideration of the energy-containing 
eddies, whose statistics are not universal. For example, the transfer functions rj + (k\k m ) and 
F + {k\k m ) obtained and plotted by Chasnov jl8| (see his Fig.l for EDQNM results and Figs. 3, 5 
for LES) cannot be transferred from one simulation to another without regard of the energy-scale 
statistics. Even such a trivial change as adding a homogeneous (zero-wavenumber) component 
to the velocity field with independent Gaussian statistics of large variance vq will drastically 
increase the size of the "spikes" in these functions at the cutoff wavenumber k m . While it 
possible to laboriously incorporate such effects with use of the sharp Fourier filter, it is totally 
unnecessary to do so when using using a sensible "good" filter. 

On the other hand, no negative conclusions can be drawn from our work regarding the 



"spectral LES" methods, such as those of Chollet and Lesieur [30], Metais and Lesieur [31], 
and others. In these LES schemes, a Fourier-Galerkin truncation is used rather than a grid- 
discretization, the stress model is a simple fc-dependent eddy-viscosity model, and the model 
dynamics is solved by a spectral code. The functional form of the eddy- viscosity in k is generally 
taken from the results of analytical closures, such as Kraichnan's fundamental work on the TFM 
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equations |3^]. Along with Chasnov's, these LES models have a somewhat different status than 
the ones we have analyzed above, based upon smooth filtering in space and a grid-discretization. 
In the latter schemes we found that the stresses in the common models (Smagorinsky, Bardina, 
etc.) will correctly reflect the magnitude of the true stresses obtained when a "good" filter 
is applied to the individual flow fields (see Eqs. (p2|),(|37|).) What we can definitely state 
regarding the "spectral LES" technique is that its stress model via the fc-dependent eddy- 
viscosity cannot be regarded as an accurate representation — realization by realization — of the 
true stress obtained by applying the sharp spectral filter to individual velocity fields in the 
turbulent ensemble. In fact, we have seen that the latter operation gives rise to a subgrid stress 
which is dominated by distant triads in individual realizations, whose effects only cancel when 
additional (scale- or ensemble-) averaging is performed. On the contrary, the eddy-viscosity 
models from the closures incorporate only the effects of the local interactions on transfer across 
the cutoff wavenumber and do not include the nonlocal transfers discussed above, which are 
actually present in individual realizations with the sharp Fourier filter! This does not rule out 
that such spectral-LES schemes will give good results, particularly for statistical averages, as 
has been often found in practice. 

On the other hand, by examination of the expressions for damping and backscatter in 
Chasnov's work [p^| , r](k\k m ) and F(k\k m ), given by his Eqs. (22) and (23), we see that the 
diverging, convective contributions are contained in these separate terms. This arises precisely 
from the range of small q in the wavenumber integrations there, i.e. q < fco an d p < k + ko. One 
expects from our analysis of the true Navier-Stokes dynamics that such contributions indeed 
ought to appear with the sharp Fourier filter but that their net contribution to mean energy 
transfer will cancel. This was explicitly verified by Chasnov through an asymptotic analysis 
of his equations (see Eqs. (20), (21) in J18||.) Hence, his stochastic model provides a really more 
faithful account of the behavior in individual realizations with the sharp Fourier filter than the 
pure viscosity- type models, as |3C], |3l]] . However, we argue that this makes the model less useful 
not more so, because it is less economical to compute and very fundamentally restricted in its 
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universality. Even for the idealized situation of homogeneous turbulence the basic quantities 
in the model must be recomputed for every new choice of low-wavenumber driving mechanism 
which leads to altered statistics of the energy-range eddies. All of these problems are totally 
eliminated by the simple expedient of using a "good" filter as defined by our criterion Eq.(p2|). 

We should emphasize again that the verb "use" has here a subtle interpretation, since in 
many LES simulations the filter is only implicit. It is often argued that the cutoff filter should 
be used to generate resolved fields that mimic those from LES with a spectral calculation, while 
the top-hat filter is more appropriate for a finite-difference scheme. However, this correspon- 



dence between numerical method and corresponding filter type is not at all exact [33]. Hence, 
whatever the numerical scheme, one is better off to employ a stress model without "connective 
backscatter." This can always be justified on the grounds that the implicit filter is a "good" 
one for which those effects are absent according to our demonstrations above. The prohibition 
against the "bad" filters, like the Fourier cutoff, is more direct if the filter is explicit in the 
model, as for the Bardina model generally or for the dynamic model. In that case, use of the 
sharp Fourier filter will require for consistency that the convective effects be included, with all 
the problems that entails. Coupled with the difficulty of using spectral codes at all in more 
complex, inhomogeneous geometries — which is a current trend of LES research — we see the 
role of the cutoff filter in the future as rather restricted. 



4 Conclusions 

Already from our work in I we can see that the numerical LES modelling technique has a 
great deal of a priori theoretical support. LES is mathematically justified even for the "weak 
solutions" of Euler equations necessary to describe non-differentiable "multifractal" velocity- 
fields in the inertial range of scales. Furthermore, many of the most common stress models 
satisfy a priori estimates obeyed by the true stress and are able to capture some of the basic 
features of inertial-range intermittency. When a "good" choice of filter function is made, the 
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subscale stress in the filtering approach is a function solely of the small scales and may be 
plausibly modelled in a universal way. However, if a "bad" filter is used — of which the cutoff 
filter is the most prominent example — then the "subgrid stress" is infected with dependence 
upon the large-scale modes and can no longer be represented in a universal form. 

One caveat concerning the LES method is that it is only really justifiable to use when 
an extended inertial interval exists. In |34j Brasseur and Wei have made an estimate of the 
Taylor-scale Reynolds number required to produce A decades of inertial range, as 

Re A ~ 400 x 10 2A/3 . (46) 

In particular, for even a single decade, He\ ~ 400 must be achieved. If this pessimistic as- 
sessment is correct, then the LES method may only begin to have a sound theoretical basis in 
circumstances of exceedingly high Reynolds number not often encountered. Of course, there is 
nothing to prevent the application of the LES method in situations where its theoretical justi- 
fication is weak and, indeed, it may do rather well there. It would not be the first time that 
the performance of a practical numerical method exceeded the parameters of its derivation. 
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